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ABSTRACT 


ERA-40 data are analyzed to demonstrate that wave forcing at lower latitudes plays a 
crucial role in driving the tropical upwelling portion of the Brewer-Dobson circulation. It is 
shown that subtropical wave forcing is correlated with tropical upwelling on both intraseasonal 
and interannual time scales when transient waves are taken into account, and that tropical wave 
forcing exerts its influence on tropical upwelling via its body force on the zonal mean flow. 


2 



1. Introduction 


Since the pioneering works of Brewer (1949) and Dobson (1956), the stratospheric meridional 
circulation, known as the Brewer-Dobson circulation (BDC), has been recognized as an integral part 
of the global climate system (Murgatroyd and Singleton 1961; Dunkerton 1978; Rosenlof and Holton 
1993; Holton et al. 1995; Haynes 2005). Model simulations show the strengthening of the BDC in 
response to greenhouse-gas induced climate change (Rind et al. 1990, 1998; Butchart et al. 2000; 
Butchart and Scaife 2001; Sigmond et al. 2004; Butchart et al. 2006, 2010, 2011; Fomichev et al. 
2007; Li et al. 2008; Garcia and Randel 2008). 

Adiabatic cooling in the upwelling branch of the BDC leads to the “tropical cold point” while 
adiabatic warming in the downwelling branch keeps temperatures above radiative equilibrium over the 
winter polar region (Andrews et al. 1987; Yulaeva et al. 1994; Salby and Callaghan 2002; Ueyama and 
Wallace 2010). Variations of the BDC have important implications for stratospheric ozone and its 
recovery (WMO/UNEP 2010 and references therein), water vapor concentrations and other chemical 
species (Ko et al. 1985; Holton 1990; Mote et al. 1996; Randel et al. 1999, 2006; Zhou et al. 2001) and 
the age of stratospheric air (Hall and Plumb 1994; Waugh and Hall 2002; Austin and Li 2006; Austin 
et al. 2007; Oman et al. 2009; Waugh 2009; Garcia et al. 2011). The extratropical part of the BDC is 
relatively better understood than tropical upwelling, i.e., the tropical part of the BDC. Wave 
dissipation in the stratosphere drives the poleward transport and forces mean downwelling at high 
latitudes (McIntyre and Palmer 1983, 1984; Haynes and McIntyre 1987; McIntyre 1990, 1999; Haynes 
et al. 1991; Holton et al. 1995; Waugh 1996), which is capsulized as the “downward control” principle 
(Haynes et al. 1991), i.e., in the steady state, the extratropical residual mean vertical velocity on a 
given pressure level is controlled exclusively by the distribution of the wave drag above that level, i.e., 
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where w* is the residual mean vertical velocity, cp the latitude, z the log pressure height, p the 
reference density profde, a the radius of the earth, D the wave forcing, fit = a cos (p(u + Qacoscp) the 
zonal mean angular momentum, u the zonal mean zonal wind, Q the angular velocity of the earth, dz' 
the vertical projection of a segment of the m curve, subscript notation for derivatives is used, and the 
integration is along angular momentum contours. 

Since conservation of mass requires that the total upwelling mass fluxes compensate the total 
downwelling mass fluxes on any pressure level, the principle of “downward control” has been 
exploited to derive mass fluxes from the troposphere to the stratosphere in the tropics and an area- 
averaged tropical upwelling velocity at 100 hPa by calculating downward mass fluxes across the 100 
hPa level in the extratropics of both hemispheres (Holton 1990; Rosenlof and Holton 1993; Rosenlof 
1995; Randel et al. 2002, 2008). Before the seminal work of Haynes et al. (1991), the residual 
velocities were only calculated from thermodynamic balance (Murgatroyd and Singleton 1961; 
Solomon et al. 1986; Gille et al. 1987; Rosenfield 1987; Shine 1989). Later on, upwelling velocities 
could also be estimated via aircraft (Boering et al. 1996) and satellite (Mote et al. 1996; Niwano et al. 
2003) measurements of trace constituents. 

However, consensus is not reached on how wave forcing drives tropical upwelling in the lower 
stratosphere. Some studies (Yulaeva et al. 1994; Holton et al. 1995; Holton 2004) argued that tropical 
upwelling is remotely driven by high-latitude wave forcing, which was termed as the ‘extratropical 
suction pump’ (Holton et al. 1995). Meanwhile, Rosenlof (1995) pointed out that tropical stratospheric 
upwelling is directly related with the wave forcing at the ‘turnaround latitudes’, i.e., the latitudes 
where the residual mean vertical velocity switches from upward to downward. The importance of 
subtropical wave forcing for tropical upwelling has been substantiated by subsequent works (Plumb 
and Eluszkiewicz 1999; Butchart and Scaife 2001; Semeniuk and Shepherd 2001; Scott 2002; 
Eichelberger and Hartmann 2005; Zhou et al. 2006; Geller et al. 2008; Randel et al. 2008; Garcia and 
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Randel 2008; McLandress and Shepherd 2009; Garny et al. 2011; Chen and Sun 2011). Other studies 
(Boehm and Lee 2003; Kerr-Munslow and Norton 2006; Norton 2006; Ryu and Lee 2010; Gamy et al. 
2011) indicate that equatorial wave forcing also plays an important role in driving tropical upwelling 
in the lower stratosphere. 

In particular, Ueyama and Wallace (2010) found that tropical upwelling is not significantly 
correlated with lower-latitude wave forcing on either intraseasonal or interannual time scales, while it 
is significantly correlated with high-latitude wave forcing on those time scales, thus casting some 
doubt on the role of lower-latitude wave forcing in driving tropical upwelling that had been established 
by the previously mentioned works. (This was mentioned in WM O/UNEP, 2011, as an issue that 
needed to be resolved by further research. See pages 4.15-4.16). Ueyama and Wallace (2010) used 
both ERA-40 monthly mean and NCEP/NCAR (Kalnay et al. 1996) pentad-mean data in their work. 
Given that a number of authors claim that both transient and stationary waves are important in forcing 
tropical upwelling, we will see if Ueyama and Wallace’s (2010) results still hold when six-hourly 
ERA-40 information is used. Also, Ueyama and Wallace (2010) found that the ‘turnaround latitude’ 
between tropical upwelling and extratropical downwelling was ~45 °, whereas previous authors found 
the ‘turnaround latitude’ to be closer to 30°. Given that Ueyama and Wallace (2010) used 
MSU/AMSU data, and that the weighting function for these is quite broad in altitude, we will examine 
the altitude structure of this ‘turnaround latitude’. 

The remainder of this paper is organized as follows. Section 2 describes the methodology and data 
used in this study. Section 3 demonstrates that both subtropical and tropical wave forcing are 
important in determining tropical upwelling. Conclusions and discussion are presented in section 4. 

2, Methodology and data 

Ueyama and Wallace (2010) used both monthly mean data and pentad-mean data to calculate 
quadratic quantities related to the Eliassen-Pahn flux and were unable to find any significant 
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correlation between tropical upwelling and lower-latitude wave forcing. They demonstrated that the 
meridional eddy heat flux computed from monthly mean ERA-40 data is similar to that calculated 
from pentad-mean NCEP-NCAR data and concluded that both datasets are reliable and give consistent 
results. However, in previous studies showing the importance of lower-latitude waves for tropical 
upwelling, either daily data (Boehm and Lee 2003; Randel 2008) or six hourly data (Kerr-Munslow 
and Norton 2006) were used to evaluate those quadratic quantities and compile monthly means. Using 
the ‘downward control’ principle and ERA-Interim data (Dee et al. 2011) from 1989 to 2008, Garny et 
al. (2011) recently showed that transient waves and stationary waves contribute comparably to tropical 
upwelling averaged between 20°S and 20°N at 100 hPa (see their Fig. 4). Randel et al. (2008) also 
pointed out that transient eddies in both hemispheres over 20°-50° contribute to tropical upwelling. 
Thus, the above-mentioned contradiction might result from missing transient waves in the monthly 
mean and pentad-mean data used by Ueyama and Wallace (2010). 

In this study, European Centre for Medium-Range Weather Forecast (ECMWF) Re -Analysis (ERA- 
40; Uppala et al. 2005) six hourly data between 1979 and 2002 are employed due to their length and 
availability up to the 1-hPa level. The horizontal resolution of the data is 2.5° latitude by 2.5° 
longitude. Data before 1979 are not used because stratospheric satellite data were not included in the 
analysis for that period. 

Ueyama and Wallace (2010) showed that temporal variations in radiative forcing, except from 
volcanic eruptions, are strictly seasonal and project almost exclusively onto the equatorial asymmetric 
component, and detennined that the temperature seesaw between the tropics and extratropics projects 
almost exclusively onto the equatorially symmetric component. In this paper, our approach follows 
closely that of Ueyama and Wallace (2010). That is to say, empirical orthogonal function (EOF) 
analysis of the equatorially symmetric component of zonally averaged monthly mean temperature in 
the lower stratosphere is used to isolate the dominant modes of variability in the BDC. 
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Since ERA-40 data are on pressure levels, the wave forcing D in (1), the Eliassen-Pahn flux F , and 
its divergence are formulated in pressure coordinates on a sphere as in Andrews et al. (1983), 


D = V-F, 

a cos cp 

F = {F (p) ,F (/,) } = acos <p{-u' v' + y/ u ,-u'a>'-y/[(a cos <p)\u cos (p\- 

- - d r)F^ p - 

V-F = ( acoscpY 1 — ( F {tp) cos^fH 

d(p dp 

_ x-T fiT 

where y/ = v'O' 1 0 = -vT/( ) . 

p dp 


( 2 ) 

(3) 

(4) 

(5) 


Here p denotes pressure, (u,v,a>) "velocity" in (longitude, latitude, pressure) coordinates, 0 potential 
temperature, T temperature, / the Coriolis parameter, and k the ratio of the gas constant to the 
specific heat at constant pressure. Primes denote departures from the zonal means, which are 
represented by overbars, and all other notations are as commonly used (Andrews et al. 1987). 

In this paper, wave forcing D in (2) is calculated via six hourly data and collected into its monthly 
mean value for analysis so that both stationary and transient waves are properly accounted for. Randel 
et al. (2004) have shown that differences exist among various stratospheric data sets. Discrepancies 
are shown in the primary observed variables (i.e., winds and temperatures), and these differences 
become larger for quadratics of these primary quantities (i.e., heat and momentum fluxes). The 
Eliassen-Pahn fluxes involve these quadratic terms, and the wave forcing involves differences in 
derivatives of these quadratic fluxes (see equations (2)-(5)). Furthermore, these differences between 
quadratic quantities involve a great deal of cancellation between terms, so the wave forcing is very 
difficult to determine accurately from reanalysis data. Geller et al. (2008) showed that close to the 
equator, small divergences of the Eliassen-Pahn flux become increasingly important. Also, F (p) 
associated with vertical eddy fluxes u'oJ in the Tropics has significant contribution from gravity 
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waves, some of which are likely not well resolved in the ERA-40 reanalysis (e. g., Bergman and Salby 
1994; Giorgetta et ah, 2002, Watanabe et al. 2008). Thus, the accuracy of the calculated EP flux 
divergences is further limited in the deep Tropics. 

In order to test the statistical significance of the correlation coefficients r(x,y) between the BDC 
and wave forcings, we use the methodology adopted by Oort and Yienger (1996) as follows. We use 
the null hypothesis that the population value p(x,y) = 0. The confidence limits were calculated via 
Fisher’s Z transformation (Spiegel et al. 2008): 



( 6 ) 



Ueyama and Wallace (2010) demonstrated that the temperature seesaw between the Tropics and 
extratropics is almost entirely captured by the even (equatorially symmetric) component of the 
meridional temperature distribution, defined as the mean of zonally averaged temperatures at 
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corresponding latitudes in the Northern and Southern Hemisphere, {T m +T SH )/2. In order to exclude 

the contribution of volcanic eruptions and the long-term cooling trend to the interannual variability 
(Ueyama and Wallace 2010), the global mean temperatures are subtracted from the meridional 
temperature distributions. The resultant meridional temperature distributions are further processed, 
with their climatological annual mean distribution (i.e., averaged over 264 months between January 
1980 and December 2001) further removed to highlight the temperature seesaw structure associated 
with the variability of wave forcing. A standard empirical orthogonal function (EOF) analysis is 
performed on these anomaly distributions weighted by the square root of the cosine of the midpoint 
latitude for each zonal belt so as to analyze area-weighted variances. Henceforth, we refer to the 
weighted anomaly distributions as T* . Note that the data in 1979 and 2002 are only used to construct 
the low-pass filtered fields for studying the interannual variability. 

For T* at the 70-hPa level, the three largest eigenvalues account for 90.1%, 5.1%, and 3.1% of the 
total month-to-month variance. Fig. 1 shows the first leading EOF and time series of its coefficient. 
The first leading EOF pattern clearly represents the temperature seesaw associated with variations of 
the BDC, i.e., the BDC in the lower stratosphere is strongest in the boreal winter and weakest in the 
boreal summer, which is consistent with the results of Yulaeva et al. (1994) among others. The 
‘turnaround latitude’ in Fig. 1(a) is approximately 3 7. 5 0 , which is comparable to the value derived 
from UKMO analysis and the results from various GCMs (Butchart et al. 2006, see their Fig. 1 and 
Fig. 2), whereas Ueyama and Wallace (2010) found the crossing latitude to be -45° using 
MSU/AMSU data. Combining EOF 1 and the corresponding first PC (principal component) in Fig. 1, 
we can infer that equatorial temperatures are above the climatological annual mean by about 2 - 4° C 
during boreal summers with maxima in Jul/Aug, and below the climatological annual mean by about 
2-4°C during boreal winters with minima in Jan/Feb, indicating that the amplitudes of the annual 
cycles in equatorial temperatures at the 70-hPa are larger than those of the equatorial temperatures 
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derived from the MSU/AMSU data (~ 2°C) by Ueyama and Wallace (2010). Fig. 1 also shows that 
the amplitude of the temperature seesaw exhibits strong interannual variability. For instance, the 
minimum equatorial temperature that occurred in February 2000 is slightly less than 2 °C below the 
climatological annual mean while that occurred in January 1998 is nearly 5 °C below the 
climatological annual mean. In addition, PC 1 shows appreciable intraseasonal variability during 
northern winters due to the stratospheric vacillations (Holton and Mass 1976; Yoden 1987) associated 
with the stronger wave forcings in those seasons. 

Similar EOF analysis is further performed on T* at the 150-hPa and 30-hPa levels to understand 
why the ‘turnaround latitude’ and the amplitudes of the annual cycles shown in Fig. 1 differ from 
those shown in Ueyama and Wallace (2001) where the EOF analysis is perfonned on the MSU/AMSU 
deep-layer atmospheric temperatures rather than the temperatures at any individual pressure level. The 
weighting function of the lower-stratospheric channel is concentrated mainly in the 30-150 hPa layer, 
as shown in Fueglistaler et al. (2011) (refer to their Fig. 3). For T* at the 150-hPa level, the three 
largest eigenvalues account for 86.9%, 7.1%, and 3.7% of the total month-to-month variance. Fig. 2 
shows the first leading EOF and time series of its coefficient at the 150-hPa level. The first leading 
EOF pattern also shows the temperature seesaw structure at the 150-hPa level. However, the 
‘turnaround latitude’ at that level is about 51° while the amplitudes of the annual cycles can be 
inferred as only about 1 U C. For T at the 30-hPa level, the three largest eigenvalues account for 
75.5%, 17.1%, and 5.8% of the total month-to-month variance. Fig. 3 shows the first leading EOF and 
time series of its coefficient at the 30-hPa level. The first leading EOF pattern at the 30-hPa level again 
represents the temperature seesaw between low latitudes and high latitudes. Both the ‘turnaround 
latitude’ and the amplitudes of the annual cycles are comparable to those in Ueyama and Wallace 
( 2010 ). 
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Reed and Vlcek (1969) showed that the annual temperature variation in the lower tropical 
stratosphere has its maximum amplitude at the 80-hPa level. It follows that the amplitudes of the 
annual cycles in the equatorial temperatures at the 70-hPa level shown in Fig. 1 are larger than those at 
the 150-hPa and 30-hPa levels shown in Fig.2 and Fig. 3 respectively. Consequently, the annual 
temperature variation at the 70-hPa level has larger amplitude than that for the 30-150 hPa layer 
shown by Ueyama and Wallace (2010). 

Previous studies (Butchart et al. 2006; McLandress and Shepherd 2009; Simpson et al. 2011; Chen 
and Sun 2011) show that the ‘turnaround latitude’ shifts poleward with increasing height below 30- 
hPa in the lower stratosphere. We found that the ‘turnaround latitude’ at the 70-hPa level shown in 
Fig. 1 is located equatorward of that at the 50-hPa level (42.5° latitude, figure not shown), which is in 
turn located equatorward of that at the 30-hPa level shown in Fig. 3. The ‘turnaround latitudes’ shown 
in Fig. 1, Fig.2, and Fig. 3 are consistent with the result of Ueyama and Wallace (2010). 

For T* at the 30-hPa level, the first eigenvalue only accounts for 75.5% of the total month-to-month 
variance while the first eigenvalues at the 70-hPa and 150-Pa levels account for 90.1% and 86.9% of 
the variances respectively, indicating that although the annual variation is still a dominant dynamical 
mode of the BDC at the 30-hPa level, other dynamical modes also play a role at this level. Since the 
second and third eigenvalues at the 30-hPa level account for 17.1%, and 5.8% respectively, the second 
EOF at the 30-hPa level is physically relevant according to North’s rule of thumb (North et al. 1982). 
Fig. 4 shows the second leading EOF and time series of its coefficient for T* at the 30-hPa level. EOF 
2 in Fig. 4(a) shows a tripolar structure with a negative-positive-negative structure with increasing 
latitude. The lower-latitude seesaw in Fig. 4(a) depicts the anticorrelation between equatorial and 
subtropical/extratropical temperatures at the 30-hPa level. This latitudinal temperature gradient is 
associated with the vertical shear of the zonal mean zonal wind at that level as shown in Fig. 4(b), due 
to the thermal wind relationship. PC 2 in Fig. 4(b) shows that the second dynamical mode is 
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dominated by the quasi-biennial oscillation (QBO). Combining EOF 2 and the corresponding PC 2 in 
Fig. 4, we can infer that the amplitudes of equatorial temperature QBO are about 4°C, which is 
comparable to those derived from the UKMO analyses by Randel et al. (1999, see their Fig. 2). Since 
the QBO is most prominent between 5-hPa and 40-hPa (Baldwin et al. 2001), EOF 1, which is 
dominated by annual variations, only accounts for 75.5% of the total month-to-month variance at the 
30-hPa level while it accounts for -90% of the total month-to-month variance at the 70-hPa and 150- 
hPa levels. This low-latitude temperature seesaw was first deduced by Reed (1964, see his Figure 8). It 
is associated with the QBO-induced meridional circulation envisaged by Reed (1964) and simulated 
by Plumb and Bell (1982). The high-latitude seesaw in Fig. 4(a) depicts the anticorrelation between 
middle-latitude and polar temperatures at the 30-hPa level, and is largely ascribed to the finding by 
Labitzke (1982) that the 30-hPa polar vortex poleward of 55°iV is much colder when the 50-hPa 
equatorial winds (see Fig. 4(b)) are from the west as per the Holton-Tan (1980, 1982) mechanism. The 
tripolar structure shown in Fig. 4(a) confirms the finding by Ruzmaikin et al. (2005) that the 
mechanism of the QBO influence on the extratropics operates mostly through the modulation of the 
northern annular mode (Thompson and Wallace 2000), However, the QBO influence on the winter 
polar vortex is an unsettled and complicated issue (Wallace and Chang 1982; Dunkerton and Baldwin 
1991; Gray et al. 2004; Naoe and Shibata 2010), and further discussion on this is beyond the scope of 
this paper. 

In summary, in this section we have shown that the leading mode of variability in the BDC between 
the 1 50-hPa and 30-hPa levels is dominated by the annual cycle with extrema in Jan/Feb and Jul/Aug, 
and that the ‘turnaround latitude’ of the BDC varies with altitude. In addition, the BDC exhibits both 
intraseasonal and interannual variability. 

4. Relationship of the BDC to wave forcing 
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Many studies have established that wave forcing at lower latitudes plays a crucial role in driving the 
tropical upwelling portion of the Brewer-Dobson circulation (Plumb and Eluszkiewicz 1999; 
Semeniuk and Shepherd 2001; Scott 2002; Eichelberger and Hartmann 2005; Zhou et al. 2006; Geller 
et al. 2008; Randel et al. 2008; Garcia and Randel 2008; Garny et al. 2011; Chen and Sun 2011). 
Ueyama and Wallace (2010) argued that if the lower-latitude wave forcing does play an important role 
in the annual cycle in the BDC as per those previous studies, then the nonseasonal variability in the EP 
fluxes might force some of the nonseasonal variability in the BDC as well. Using pentad NCEP- 
NCAR data and monthly mean ERA-40 data, Ueyama and Wallace (2010) found that the nonseasonal 
variability in the BDC is not significantly correlated with that in lower-latitude wave forcing. In this 
section, we show that nonseasonal variabilities in the BDC and the lower wave forcing are indeed 
significantly correlated if six-hourly ERA-40 data are employed, thus properly taking transient waves 
into account. 


It is pointed out in the section 2 that the accuracy of the calculated EP flux divergences is very 
limited in the deep Tropics. In addition, the ‘downward control’ principle is not applicable in the deep 
Tropics due to the fact that the gradient of angular momentum, i.e., m in Eqn. (1), tends to vanish 

there (Haynes et al. 1991). Thus, we only present how the BDC is correlated to wave forcings between 
15° and 85° latitude where Eqn. (1) can be approximated (Haynes et al. 1991) by 


w\<p,z) = - 


1 d cos (p 

pa cos (p dcp 2Qsin#> 


J *oo 
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Combining (2) and (10) leads to 
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where g is the acceleration of gravity and other symbols are described in previous sections. 
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Similarly to Rosenlof 1995 and Butchart et al. 2006, the 70-hPa level is selected to represent the 
lower stratosphere for purpose of investigating the BDC. The corresponding wave forcing at each 
parallel is measured by 

WF =\ P \V -F) , dp (12) 

J Pi T 

where p x =\hPa and p 2 =70hPa . 

a. Nonseasonal variabilities of the lower stratospheric BDC and wave forcing 

In order to investigate nonseasonal variabilities of the lower stratospheric BDC and wave forcing, 
the climatological mean annual cycles over the period 1980-2001 are further removed from T* at the 
70-hPa level. The corresponding wave forcing (WF expressed in eqn. (12)), is treated in the same 
manner, so that associations between the deseasonalized T* at the 70-hPa level and the deseasonalized 
WF can be investigated. 

An EOF analysis was performed on the deseasoned T * at the 70-hPa level. The largest eigenvalue 
accounts for 57.5% of the total month-to-month variance. Fig. 5 shows the first leading EOF and time 
series of its coefficient. The first leading EOF pattern still represents the temperature seesaw 
associated with nonseasonal variations of the BDC in the lower stratosphere with its ‘turnaround 
latitude’ located near 37.5°, which is ~ 10° equatorward of that derived by Ueyama and Wallace 
(2010). As indicated earlier, their value represents the mean crossing latitude over the 30-150 hPa 
layer. PC 1 and EOF 1 in Fig. 5 indicate that the amplitudes of the nonseasonal variations in equatorial 
temperatures at the 70-hPa are very irregular with the maximum amplitude being ~ 3°C . 

Fig. 6 shows the latitudinal distribution of the contemporary correlation coefficients between the PC 
1 of the deseasonalized T* at the 70-hPa level and the deseasonalized WF expressed in (12). Fig. 6(a) 
illustrates that the wave forcing between 17.5° and 65° latitude coherently, significantly and 
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negatively correlates with PC 1 at the 95% confidence level when six-hourly ERA-40 data are used. 
However, Fig. 6(b) shows that wave forcing in the subtropics no longer significantly correlates with 
the PC 1 at the 95% confidence level when monthly mean ERA -40 data are used even though the 
wave forcing between 42.5° and 67.5° latitude still significantly correlates with the PC 1 at the 95% 
confidence level. 

Fig. 6 shows that the PCI are negatively and significantly correlated with wave forcings between 
42.5° and 65° latitude no matter what kind of dataset is used, which confirms the importance of 
planetary- wave breaking in the stratospheric ‘surf zone’ (McIntyre and Palmer 1983, 1984) to the 
BDC. Fig. 6(b) is also largely consistent with the results of Ueyama and Wallace (2010), namely, 
tropical upwelling in the BDC is significantly correlated with the high-latitude wave forcing. Contrary 
to Ueyama and Wallace (2010), however, Fig. 6 (a) clearly demonstrates that subtropical wave forcing 
plays an important role in driving the tropical upwelling portion of the Brewer-Dobson circulation as 
shown by Plumb and Eluszkiewicz (1999), Zhou et al. (2006), Geller et al. (2008) and Chen and Sun 
( 2011 ). 

Fig. 6 (a) also shows that the PCI are positively and significantly correlated with wave forcings 
between 72.5° and 82.5° latitude. However, Geller et al. (2008) pointed out that tropical upwelling in 
the BDC is not sensitive to the wave forcing in the polar region, thus we do not further look into the 
role of polar wave forcing. 

b. Interannual variabilities of the lower stratospheric BDC and wave forcing 

In order to investigate interannual variabilities of the lower stratospheric BDC and wave forcing, 
we further apply successively centered five- and three-month running mean smoothing operators to the 
deseasonalized T* at the 70-hPa level and the deseasonalized wave forcing to form the low-pass 
filtered version (Ueyama and Wallace 2010). Fig. 7 (a) shows EOF 1 of the low-pass filtered and 
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deseasonalized T * , which accounts for 58.2% of the total month-to-month variance. In the low-pass 

filtered pattern, the crossing latitude is located near 22.5°, which is - 15° equatorward of the mean 
crossing latitude over the 30-150 hPa layer derived by Ueyama and Wallace (2010). PC 1 and EOF 1 
in Fig. 7 indicate that the amplitude of the interannual variations in equatorial temperatures at the 70- 
hPa is about 2°C. PC 1 apparently exhibits a QBO variation. Fig. 7 (b) shows that PC 1 is negatively 
correlated with the vertical shear of the zonal mean equatorial zonal wind at the 70-hPa level 
represented by u 50hPa - u 10hPa . The contemporary correlation coefficient is -0.57. The degree of 

freedom is 16. Thus, the two-tailed test verifies that the vertical shear and PC 1 are significantly 
correlated at the 95% confidence level. 

Fig. 8 shows the latitudinal distribution of the contemporary correlation coefficients between the PC 
1 of the low-pass filtered and deseasoned T* at the 70-hPa level and the low-pass filtered and 
deseasoned WF expressed in (12). Fig. 8 (a) demonstrates that wave forcing between 20° and 62.5° 
latitude coherently, significantly and negatively correlates with PC 1 at 95% confidence level when 
six-hourly ERA-40 data are used. However, Fig. 8(b) shows that wave forcing in the subtropics no 
longer significantly correlates with the PC 1 at 95% confidence level when monthly mean ERA-40 
data are used even though wave forcing around 45° and 65° latitude respectively still significantly 
correlates with the PC 1 at the 95% confidence level. Thus, we conclude that subtropical wave forcing 
does play an important role in driving the tropical upwelling portion of the Brewer-Dobson circulation 
on interannual time scale when the contributions from transient waves are considered. 

Salby (2011) pointed out that the variation of wave forcing alone cannot reproduce the observed 
response of stratospheric temperature and ozone in his model, and that nearly all of the interannual 
variance of stratospheric temperature and ozone over the Northern Hemisphere can be accounted for if 
both the variation of wave forcing and the QBO of equatorial wind are included. 
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c. Intraseasonal variabilities of the lower stratospheric BDC and wave forcing 

In order to investigate intraseasonal variabilities of the lower stratospheric BDC and wave forcing, 
we subtract the low-pass fdtered and deseasonalized T* at the 70-hPa level and the corresponding 
wave forcing from the deseasonalized T* and wave forcing respectively. ). Fig. 9 (a) shows EOF 1 of 
the high-pass fdtered and deseasonalized T * , which accounts for 76.0% of the total month-to-month 
variance. In the high-pass filtered pattern, the crossing latitude is located near 47.5° , which is roughly 
the same as the mean crossing latitude, i.e., 50°, over the 30-150 hPa layer derived by Ueyama and 
Wallace (2010). In addition, the pattern is also characterized by remarkably uniform amplitudes 
equatorward of the crossing latitude and much larger temperature anomalies of opposite sign in the 
polar region as pointed out by Ueyama and Wallace (2010). PC 1 and EOF 1 in Fig. 9 indicate that 
that the amplitudes of the intraseasonal variations in equatorial temperatures at the 70-hPa are very 
irregular with the maximum amplitude being ~ 1°C . 

Fig. 10 shows the latitudinal distribution of the contemporary correlation coefficients between the 
PC 1 of the high-pass filtered and deseasonalized T* at the 70-hPa level and the high-pass filtered and 
deseasonalized WF expressed in eqn. (12). Fig. 10 (a) demonstrates that the wave forcing between 
27.5° and 62.5° latitude coherently, significantly and negatively correlates with PC 1 at the 95% 
confidence level when six-hourly ERA-40 data are used. Fig. 10 (b) shows that only wave forcing 
between 55° and 65° latitude is significantly correlated with PC 1 at the 95% confidence level when 
monthly mean ERA-40 data are used. 

Equations (1) and (11) are only valid in the steady state (Haynes et al. 1991). Rosenlof and Holton 
(1993) indicated that application of the downward control principle appears to only be valid for the 
solstice seasons. They pointed out that the assumption of steady state is violated in the equinox 
seasons, leading to an inaccurate estimate of the residual circulation stream function during those 
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seasons if the downward control principle is employed. Since wave forcing drives the mean meridional 
circulation and accelerates the zonal mean zonal wind on intraseasonal time scales, the downward 
control principle (1) can be modified as follows (Randel et al. 2002, 2008). 


w\cp,z) 
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The modified version of Eqn. (11) can be written as 
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is referred to as the adjusted wave forcing. 

Fig. 1 1 shows the latitudinal distribution of the contemporary correlation coefficients between the 
PC 1 of the high-pass filtered and deseasonalized T* at the 70-hPa level and the high-pass filtered and 
deseasonalized adjusted wave forcing expressed by eqn. (15). In comparison with Fig. 10 (a), , Fig. 1 1 
(a) illustrates that wave forcing between 22.5° and 65° latitude coherently, significantly and more 
negatively correlates with PC 1 at the 95% confidence level when six-hourly ERA-40 data are used. In 
comparison with Fig. 10 (b), Fig. 1 1 (b) shows that the wave forcing between 52.5° and 67. 5° latitude 
is significantly and more negatively correlated with PC 1 at the 95% confidence level when monthly 
mean ERA-40 data are used. Thus, we conclude that subtropical wave forcing indeed plays an 
important role in driving the tropical upwelling portion of the Brewer-Dobson circulation on 
intraseasonal time scales, and that the effects of transient eddies are important and must be considered. 

4. Discussion and conclusions 
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The recent paper by Ueyama and Wallace (2010) found that tropical upwelling into the lower 
stratosphere is mainly under the control of high-latitude wave driving. This is a surprising result in 
light of previous theoretical works (e. g., Plumb and Eluszkiewicz 1999; Semeniuk and Shepherd 
2001; Zhou et al. 2006) that showed that the near-equatorial and subtropical wave driving were crucial 
in maintaining annually averaged tropical upwelling. This has motivated us to essentially redo the 
Ueyama and Wallace (2010) analysis using data sets with more frequent time sampling to see how this 
might alter their conclusions. 

The main results of our analyses are as follows: 

1) At the 70-hPa level, the temperature seesaw associated with variations of the BDC accounts for 
90.1% of the month-to-month variance of the temporal variations about the annual-mean 
latitudinal distribution. This dynamic mode is dominated by the annual cycle, i.e., the equator-to 
subpolar temperature contrast is largest in Jan/Feb and smallest in Jul/Aug. 

2) The ‘turnaround latitude’ associated with the annual cycle of the BDC shifts poleward with 
increasing height below 30-hPa in the lower stratosphere. At the 70-hPa level, it is approximately 
equal to 37.5° , which is comparable to the value derived from UKMO analysis and the results 
from various GCMs. 

3) At the 70-hPa level, the temperature seesaw associated with nonseasonal variations of the BDC 
accounts for 57.5% of the month-to-month variance of the temporal variations about the 
climatological mean annual cycle of latitudinal temperature distribution. The ‘turnaround latitude’ 
associated with nonseasonal variations of the BDC is the same as that associated with the annual 
cycle of the BDC. The nonseasonal variability of subtropical wave forcing is coherently and 
significantly correlated with that of the BDC in the lower stratosphere. 
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4) At the 70-hPa level, the temperature seesaw associated with interannual variations of the BDC 

accounts for 58.2% of the month-to-month variance of the temporal variations in the low-pass 
filtered and deseasonalized latitudinal temperature distribution. The ‘turnaround latitude’ 

associated with interannual variations of the BDC is located near 22. 5° . The interannual variability 
of subtropical wave forcing is indeed coherently and significantly correlated with that of the BDC 
in the lower stratosphere. 

5) At the 70-hPa level, the temperature seesaw associated with interannual variations of the BDC 
exhibits a QBO variation. This is consistent with the fact that its interannual variability is 
significantly correlated with the vertical shear of the equatorial mean zonal wind at that level. 

6) At the 70-hPa level, the temperature seesaw associated with intraseasonal variations of the BDC 

accounts for 76.0% of the month-to-month variance of the temporal variations in the high-pass 
filtered and deseasonalized latitudinal temperature distribution. The ‘turnaround latitude’ 

associated with intraseasonal variations of the BDC is located near 47.5° . The intraseasonal 
variability of subtropical wave forcing is physically and significantly correlated with that of the 
BDC in the lower stratosphere. 

7) Effects of transient eddies are important and must be considered. Very different conclusions are 
reached regarding latitudinal ranges over which the wave forcing (e.g., Eliassen-Palm flux 
convergence) contributes to the forcing of the tropical upwelling part of the Brewer-Dobson 
circulation when six -hourly data are used instead of monthly mean and pentad data that were used 
by Ueyama and Wallace (2010) since these infrequent data samplings do not properly resolve 
transient wave forcing. 
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8) Subtropical wave forcing plays a crucial role in driving the tropical upwelling portion of the 
Brewer-Dobson circulation on various time scales when the effects of transient waves are properly 
considered. 


We confirm Ueyama and Wallace’s (2010) finding that wave forcing in the stratospheric ‘surf zone’ 
plays an important role in driving tropical upwelling in the BDC on various time scales. However, 
since tropical upwelling in the BDC is coherently and significantly correlated with wave forcings in 
both the subtropical and the stratospheric ‘surf zone’ regions, we conclude that wave forcings in those 
regions act synergistically to drive tropical upwelling in the BDC. Although wave forcing near the 
‘turnaround latitudes’ can largely control the strength of the BDC (Rosenlof 1995; Butchart and Scaife 
2001; Fomichev et al. 2007; Shepherd and McLandress 2011), the latitudinal distribution of the 
residual mean vertical velocity, thus the ‘turnaround latitude’ per se, is determined jointly by both low- 
latitude and higher-latitude wave forcings. 

Butchart and Scaife (2001) showed that increased tropical upwelling in response to increased 
greenhouse-gas concentration induced climate change results from strengthened wave forcing in the 
subtropical lower stratosphere, which in turn is largely caused by the increased penetration of resolved 
waves into that region due to the amplified zonal winds there. Fomichev et al. (2007) also showed that 
the enhanced tropical upwelling of the BDC is linked to increases in subtropical wave forcing rather 
than changes in wave forcing at higher latitudes. 

Although climate models project the enhancement of tropical upwelling as a result of climate 
change (Butchart et al. 2006, 2010), consensus has not yet been reached on the exact manner in which 
tropical upwelling in the lower stratosphere is increased by wave forcing. Rind et al. (1990, 1998) 
argued that the increases in resolved wave forcing arise from increased generation of planetary waves 
due to decreases in large-scale static stability, and from improved propagation conditions associated 
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with changes in refractive index due to strengthened upper flank of the subtropical jets. Rind et al. 
(1990) also pointed out that increased gravity wave forcing results from greater wave generation in 
response to climate changes, and from higher breaking altitudes of gravity waves due to strengthened 
upper fla nk of the subtropical jets, leading to more gravity wave momentum fluxes transferred into the 
stratosphere, which is confirmed by Li et al. (2008) and McLandress and Shepherd (2009). However, 
Butchart and Scaife (2001) and McLandress and Shepherd (2009) showed a strengthened BDC 
without any discernible change in refractive index. 

Plumb (2002) depicted a paradigm where the BDC is driven by both planetary waves and synoptic 
waves. McLandress and Shepherd (2009) showed that together with parameterized orographic gravity 
wave forcing, both planetary waves and synoptic waves contribute to the enhancement of tropical 
upwelling. Shepherd and McLandress (2011) revealed that transient waves account for 80% of the 
contribution of resolved waves to the strengthened BDC in response to increased greenhouse-gas 
concentration induced climate change with roughly equal contributions from planetary waves and 
synoptic waves. Shepherd and McLandress (2011) further demonstrated that strengthened upper flank 
of the subtropical jets leads to the uplifting of the critical-layer associated with Rossby-wave breaking 
(Randel and Held 1991), which in turn makes for more Rossby wave activity penetrating into the 
subtropical lower stratosphere, giving rise to strengthened transient wave forcing there. Thus, critical- 
layer control of wave breaking provides a unified theoretical framework within which gravity waves 
and transient waves of both planetary and synoptic scale give rise to enhanced subtropical wave 
forcing in the lower stratosphere in response to increased greenhouse-gas concentration induced 
climate change. 

Although reanalysis data in the deep tropics are not reliable to evaluate the wave forcing there, we 
have found that the interannual variability of the BDC is significantly correlated with the vertical shear 
of the equatorial mean zonal wind. Wallace and Holton (1968) and Holton (1968) conclusively 
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demonstrated that tropical wave forcing is required to account for the QBO of the zonal wind in the 
tropical stratosphere. To say the least, tropical wave forcing plays an indirect role in driving the 
tropical upwelling portion of the Brewer-Dobson circulation through its action on the zonal mean 
flow. Kerr-Munslow and Norton (2006) and Norton (2006) showed that tropical stationary waves 
drive the annual cycle in tropical tropopause temperatures while Gamy et al. (2011) showed that 
tropical stationary waves play a leading role in the acceleration of upwelling velocities in the deep 
tropics in response to increases in greenhouse gases. 

Thus, we conclude that both low-latitude and high latitude wave forcing are important in driving 
tropical upwelling in the Brewer-Dobson circulation. In particular, the relative strengths of these 
forcings are important in determining the latitudinal distribution of the tropical upwelling. 
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Figure captions 


Fig. 1. The first leading EOF of T* at the 70-hPa level from 1980 to 2001 (a) and the corresponding 
first principal component (b). The horizontal scale in (a) is such that the spacing between latitudes is 
proportional to the area of the earth’s surface between them, i.e., is linear in the sine of the latitude 
(Gill 1982). 

Fig. 2. Same as Fig. 1 , but for T* at the 150-hPa level from 1980 to 2001. 

Fig. 3. Same as Fig. 1 , but for T * at the 30-hPa level from 1980 to 2001. 

Fig. 4. (a) The second leading EOF of T at the 30-hPa level from 1980 to 2001. The horizontal scale 
is the same as in Fig. 1(a). (b) The corresponding second principal component (black line), the vertical 
shear of the zonal mean equatorial zonal wind at the 30-hPa level represented by u 20hPa - u V)hPa (red 
line), and the zonal mean equatorial zonal wind at the 50-hPa level (blue line). 

Fig. 5. The first leading EOF of the deseasonalized T* at the 70-hPa level from 1980 to 2001 (a) and 
the corresponding first principal component (b). The horizontal scale in (a) is the same as in Fig. 1(a). 

Fig. 6. Contemporary correlation coefficients (solid line) between PC 1 of the deseasonalized T at the 
70-hPa level and the deseasonalized wave forcings from 1980 to 2001 derived from the six-hourly 
ERA-40 data (a) and the monthly mean ERA-40 data (b). Dotted lines denote the 95% confidence 
limits as a function of latitude. 

r i * 

Fig. 7. (a) The first leading EOF of the low-pass filtered and deseasonalized ‘ at the 70-hPa level 
from 1980 to 2001. The horizontal scale is the same as in Fig. 1(a). (b) The corresponding first 
principal component (solid line), the vertical shear of the zonal mean equatorial zonal wind at the 70- 

hPa level represented by UsohPa u iohp a (dotted line). 

Fig. 8. Same as Fig. 6, but for the low-pass filtered and deseasonalized T* at the 70-hPa level and the 
corresponding wave forcing. 
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Fig. 9. Same as Fig. 1, but for the high-pass filtered and deseasonalized T* at the 70-hPa level. 

Fig. 10. Same as Fig. 6, but for the high-pass filtered and deseasonalized T at the 70-hPa level and the 
corresponding wave forcing. 

Fig. 11. Same as Fig. 6, but for the high-pass filtered and deseasonalized T* at the 70-hPa level and the 
corresponding adjusted wave forcing. 
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(a) EOF 1 at the 70-hPa leve 




Year 

Fig. 1. The first leading EOF of T* at the 70-hPa level from 1980 to 2001 (a) and the corresponding 
first principal component (b). The horizontal scale in (a) is such that the spacing between latitudes is 
proportional to the area of the earth’s surface between them, i.e., is linear in the sine of the latitude 
(Gill 1982). 
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EOF 2 at the 30 — hPa leve 




Fig. 4. (a) The second leading EOF of T* at the 30-hPa level from 1980 to 2001. The horizontal scale 
is the same as in Fig. 1(a). (b) The corresponding second principal component (black line), the vertical 
shear of the zonal mean equatorial zonal wind at the 30-hPa level represented by u 20hPa - u 30hPa (red 
line), and the zonal mean equatorial zonal wind at the 50-hPa level (blue line). 
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(a) EOF 1 of the deseasonalized T‘ 




Fig. 5. The first leading EOF of the deseasonalized T* at the 70-hPa level from 1980 to 2001 (a) and 
the corresponding first principal component (b). The horizontal scale in (a) is the same as in Fig. 1(a). 
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(a) Correlation coefficients derived from the six-hourly data 
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(b) Correlation coefficients derived from the monthly mean data 
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Fig. 6. Contemporary correlation coefficients (solid line) between PC 1 of the deseasonalized T* at the 
70-hPa level and the deseasonalized wave forcings from 1980 to 2001 derived from the six-hourly 
ERA-40 data (a) and the monthly mean ERA-40 data (b). Dotted lines denote the 95% confidence 
limits as a function of latitude. 
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(a) EOF 1 , low-pass 




Fig. 7. (a) The first leading EOF of the low-pass filtered and deseasonalized T * at the 70-hPa level 
from 1980 to 2001. The horizontal scale is the same as in Fig. 1(a). (b) The corresponding first 
principal component (solid line), the vertical shear of the zonal mean equatorial zonal wind at the 70- 
hPa level represented by u 50hPa ~u 10hPa (dotted line). 
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